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Abstract 



We discuss the appearance at the QCD phase transition, and the subsequent 
decay, of axion walls bounded by strings in = 1 axion models. We argue 
on intuitive grounds that the main decay mechanism is into barely relativistic 
axions. We present numerical simulations of the decay process. In these 
simulations, the decay happens immediately, in a time scale of order the 
light travel time, and the average energy of the radiated axions is {tOa) — 
7ma for Va/rria ~ 500. {u)a) is found to increase approximately linearly with 
\n.{va/ma)- Extrapolation of this behaviour yields {uj a) ~ 60 771 (J in axion 
models of interest. We find that the contribution to the cosmological energy 
density of axions from wall decay is of the same order of magnitude as that 
from vacuum realignment, with however large uncertainties. The velocity 
dispersion of axions from wall decay is found to be larger, by a factor 10^ or 
so, than that of axions from vacuum realignment and string decay. We discuss 
the implications of this for the formation and evolution of axion miniclusters 
and for the direct detection of axion dark matter on Earth. Finally we discuss 
the cosmology of axion models with > 1 in which the domain wall problem 
is solved by introducing a small Upq(1) breaking interaction. We find that 
in this case the walls decay into gravitational waves. 

PACS numbers:14.80.Mz,95.30.Cq 
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I. INTRODUCTION 



The axion |]T],3 was postulated approximately twenty years ago to explain why the strong 
interactions conserve the discrete symmetries P and CP in spite of the fact that the Standard 
Model of particle interactions as a whole violates those symmetries. It is the quasi-Nambu- 
Goldstone boson associated with the spontaneous breaking of a Upq{1) symmetry which 
Peccei and Quinn postulated. At zero temperature the axion mass is given by: 

~ 6 • lO-'^eV ■ N ■ (^^^^^^^ (1-1) 

where Va is the magnitude of the vacuum expectation value that breaks Upq{1) and is a 
strictly positive integer that describes the color anomaly of Upq{1). Axion models have A^ 
degenerate vacua . Searches for the axion in high energy and nuclear physics experiments 
have only produced negative results. By combining the constraints from these experiments 
with those from astrophysics one obtains the following bound: < 10~^ eV. 

The axion owes its mass to non-perturbative QCD effects. In the very early universe, at 
temperatures high compared to the QCD scale, these effects are suppressed and the axion 
mass is negligible. The axion mass turns on when the temperature approaches the QCD 
scale and increases till it reaches the value given in Eq.(l.l) which is valid below the QCD 
scale. There is a critical time ti, defined by ma{ti)ti = 1, when the axion mass effectively 
turns on 0. The corresponding temperature Ti ~ 1 GeV. 

The implications of the existence of an axion for the history of the early universe may be 
briefly described as follows. At a temperature of order Va, a phase transition occurs in which 
the Upq{1) symmetry becomes spontaneously broken. This is called the PQ phase transition. 
At that time axion strings appear as topological defects. One must distinguish two cases: 
1) inflation occurs with reheat temperature higher than the PQ transition temperature 
(equivalently, for the purposes of this paper, inflation does not occur at all) or 2) inflation 
occurs with reheat temperature less than the PQ transition temperature. 

In case 2) the axion field gets homogenized by inflation and the axion strings are 'blown 
away'. When the axion mass turns on at ti, the axion field starts to oscillate. The amplitude 
of this oscillation is determined by how far from zero the axion field is when the axion mass 
turns on. The axion field oscillations do not dissipate into other forms of energy and hence 
contribute to the cosmological energy density today 0. Such a contribution is called of 
"vacuum realignment" . Note that the vacuum realignment contribution may be accidentally 
suppressed in case 2) because the axion field, which has been homogenized by inflation, may 
happen to lie close to zero. 

In case 1) the axion strings radiate axions from the time of the PQ transition till 
ti when the axion mass turns on. At ti each string becomes the boundary of A^ domain 
walls. If A^ = 1, the network of walls bounded by strings is unstable P,p!0| and decays away. 
All of the present paper except section VI is concerned with this process, how many axions 
are produced in the decay, what is the velocity dispersion of these axions and what are the 
implications. 

If A^ > 1 there is a domain wall problem [0] because axion domain walls end up domi- 
nating the energy density, resulting in a universe very different from the one observed today. 
There is a way to avoid the domain wall problem by introducing an interaction which slightly 
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lowers one of the vacua with respect to the others. In that case, the lowest vacuum takes 
over after some time and the domain walls disappear. There is little room in parameter 
space for that to happen but it is a logical possibihty. Section VI discusses this in detail. 

In case 1) there are three contributions to the axion cosmo logical energy density. One 
contribution |^,^,|ri|,|l^] is from axions that were radiated by axion strings before ti; let 
us call it the string decay contribution. A second contribution is from axions that were 
produced in the decay of walls bounded by strings after ti |]T2|-Q; call it the contribution 
from wall decay. A third contribution is from vacuum realignment 0. To convince oneself 
that there is a vacuum realignment contribution distinct from the other two, consider a 
region of the universe which happens to be free of strings and domain walls. In such a 
region the axion field is generally different from zero, even though no strings or walls are 
present. After time ti, the axion field oscillates implying a contribution to the energy density 
which is neither from string decay nor wall decay. Since the axion field oscillations caused 
by vacuum realignment, string decay and wall decay are mutually incoherent, the three 
contributions to the energy density should simply be added to each other [|12 . 



That axions walls bounded by strings decay predominantly into barely relativistic axions 
was suggested in ref. [1^, where the size of the wall contribution to the axion energy density 
was estimated to be of the same order of magnitude as that from string decay and from 
vacuum realignment. This is consistent with what we find here. D. Lyth |13] also discussed 
the wall decay contribution and emphasized the uncertainties affecting it. M. Nagasawa and 
M. Kawasaki [jl4[ performed computer simulations of the decay of walls bounded by string 
and obtained {ua)/ma — 3 for the average energy of the radiated axions. Our simulations, 
presented in section IV, are similar to those of Nagasawa and Kawasaki but they are done 
on larger lattices and for a wider variety of initial conditions. We obtain {uOa)/ma — 7 for 
Va/fTia — 500. We attribute the difference between our result and that of Nagasawa and 
Kawasaki to the fact that we give angular momentum to the collapsing walls, whereas they 
did not. 

It should be emphasized that the lattice sizes which are amenable to present day com- 
puters are at any rate small compared to what one would ideally wish. Indeed the axion 
string core has size of order where A is a coupling constant (see section ||) whereas the 

axion domain wall thickness is of order m^,^ . The lattice constant must be smaller than 
for the lattice to resolve the string core. On the other hand the lattice size must be larger 
than to contain at least one wall. Hence the lattice size in units of the lattice constant 
must be of order 10 ^^° x 10 ^"° or larger if the simulations are done in 2 dimensions (see 

ma ma ° ^ 

section IV). Present day computers allow lattice sizes of order 4000 x 4000, i.e. ^ ~ 100. 

In axion models of interest ^^^^ ~ ^ify-^cW ~ The computer simulations inform us 

about the situation of interest only insofar it may be assumed that the motion and decay 
of walls bounded by strings do not vary dramatically from the case where — is large to the 
case where ^ is huge. To address this issue, we study the dependence of {u!a)/ma upon 

the and find that it increases approximately as the logarithm of that quantity over the 
range accessible to the simulations. This behaviour is understood in terms of the process by 
which the walls bounded by string decay. If it persists all the way up to ^ = 10^6, then 
{uja)/ma is of order 60 in models of interest. 

There has been disagreement in the literature about the size of the string contribution. 



3 



Many authors [0|TT1 believe it is about a factor 100 larger than the vacuum realignment 



contribution. We ourselves |g|T^ have argued that the string and vacuum realignment 
contributions are of the same order of magnitude. The two different estimates are presented 
in subsection IIIB. With regard to the wall contribution, we argue that it is again of the 
same order of magnitude as that from vacuum realignment, with however large uncertainties. 
Section 111 presents a unified treatment of all contributions. The resulting axion cosmological 
energy density is obtained in subsection HID. The requirement that it does not overdose 
the universe, implies ^ 10~^ eV. By combining this bound with the one from stellar 
evolution mentioned earlier, the axion mass is constrained to lie in the window: 10~^eV 
> m„ > 10-6 eV. 

We find that the velocity dispersion of the axions from wall decay is much larger, by a 
factor of 10^ or so, than the velocity dispersion of the axions from string decay and vacuum 
realignment. This has interesting implications for the formation and evolution of axion 
mini-clusters [rB,|T^ as discussed in section V. We find that the axions from wall decay bind 



only very loosely to axion mini-clusters and that they get readily stripped off when the mini- 
cluster falls into a galactic halo. This ensures the existence of an unclustered component of 
galactic axion dark matter. The existence of such a component is important for searches of 
axion dark matter on Earth [|TB1 because it guarantees that the signal is on at all times. 



Another interesting consequence of the larger velocity dispersion of the axions from wall 
decay is that it may conceivably be measured. If a signal is found in the cavity detector 
of dark matter axions |jl8|, the energy spectrum of these axions can be measured with 



great resolution. It has been pointed out that there are peaks in the spectrum |]19[ because 
late infall of cold dark matter onto our galaxy produces distinct flows, each one with a 
characteristic local velocity vector. These peaks are broadened by the primordial velocity 
dispersion. The minimum time required to measure the primordial velocity dispersion of 
axions from vacuum realignment and string decay is too long (~ 10 years) to be practical, 
but the minimum time required to measure that of axions from wall decay is much less. 
This is discussed in subsection lllE. 

This paper is organized as follows. In section 11, we present a toy model which captures 
all the physics we are interested in. We discuss the evolution of axion strings in the early 
universe and the formation of walls bounded by strings at the onset of the QCD phase 
transition. We give an estimate of the density of walls at that time. In section III, we discuss 
the decay of walls bounded by string and the axion population that results therefrom. We 
also discuss the axions from vacuum realignment and from string decay and estimate the 
total cosmological energy density from these mechanisms. In section IV, we present our 
numerical simulations of the motion and decay of walls bounded by strings. In section V, 
we discuss the effect of axions from wall decay upon the formation and evolution of axion 
miniclusters. In section VI, we discuss axion models in which > 1 and in which the axion 
domain wall problem is solved by slightly lowering one of the vacua with respect to the 
others. In section VII, we summarize our conclusions. 



II. THE NETWORK OF WALLS BOUNDED BY STRINGS 

The bulk of this paper is concerned with the cosmological properties of A^ = 1 axion 
models. A^ = 1 means that the axion model has a unique vacuum. The following toy model 
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incorporates the qualitative features of interest to us: 



= \d,<P,d^ct>, + \d,<P2d^ct>2 - ^(0? + 4>l- vlf + r^vact>i (2.1) 

where = 0i + i(j)2 is a complex scalar field. When 77 = 0, this model has a U{1) global 
symmetry under which — > e*"0(x). It is analogous to the Upq{1) symmetry of Peccei 
and Quinn. It is spontaneously broken by the vacuum expectation value (0) = fae*"". The 
associated Nambu-Goldstone boson is the axion. The last term in Eq. represents the 
non-perturbative QCD effects that give the axion its mass nia- We have = ^Jfj for 

When ?7 = 0, the model has global string solutions. A straight global string along the 
5-axis is the static configuration: 

m = vj{p)e'' (2.2) 

where {z,p,6) are cylindrical coordinates and /(p) is a function which goes from zero at 
p = to one at p = 00 over a distance scale of order 6 = {y/Xva}^^- f{p) may be determined 
by solving the non-linear field equations derived from Eq. ( |2.1| ). The energy per unit length 
of the global string is 

p ^ 27r pdpl I V0p = Tivl Hy^VaL) (2.3) 

where L is a long-distance cutoff. Eq. (2.3) neglects the energy per unit length, of order 
associated with the core of the string. For a network of strings with random directions, as 
would be present in the early universe, the infra-red cutoff L is of order the distance between 
strings. 

When ?7 7^ 0, the model has domain walls. If Va 3> vna we may set 0(x) = VaC.'^"'^^^ when 
restricting ourselves to low energy configurations. The corresponding effective Lagrangian 
is : 



r (a\ 




COS — 


- 1 







Its equation of motion has the well-known Sine-Gordon soliton solution : 



(2.4) 



27rn-|-4tan ^exp(ma2/) (2.5) 



where y is the coordinate perpendicular to the wall and n is any integer. Eq. ( |2.5| ) describes 
a domain wall which interpolates between neighboring vacua in the low energy effective 
theory ( ^.41 ). In the original theory ( |2.1| ), the domain wall interpolates between the unique 
vacuum back to that same vacuum going around the bottom of the Mexican hat potential 

The thickness of the wall is of order . The wall energy per unit surface is cr = 8maf^ 
in the toy model. In reality the structure of axion domain walls is more complicated than 
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in the toy model, mainly because not only the axion field but also the neutral pion field is 



spatially varying inside the wall |20]. When this is taken into account, the (zero temperature) 



wall energy/surface is found to be: 

a ^ 4.2 UmJa ^ 9 mj! (2.6) 

with fa = Va/N. For N = 1, fa and Va are the same. However, written in terms of fa, Eq. 
(2.6) is valid for 7^ 1 also. 

When rria 7^ 0, each string becomes the boundary of a domain wall. Indeed the presence 
of the string forces a{x) = ■:^a{x) to vary by 2tt when x goes around the string once and a 
domain wall is the minimum energy path that does this. Fig. 1 shows a cross-section of a 
wall bounded by a string. Two different length scales are involved: the wall thickness 
and the string core size 6 = {\/Xva)~^. Walls bounded by string are of course topologically 
unstable. 

In the early universe, the axion is essentially massless from the PQ phase transition at 
temperature of order Va to the QCD phase transition at temperature of order 1 GeV. During 
that epoch, axion strings are present as topological defects. At first the strings are stuck 
in the plasma and are stretched by the Hubble expansion. However with time the plasma 
becomes more dilute and below temperature M 



the strings move freely and decay efficiently into axions. Because this decay mechanism is 
very efficient, there is approximately one string per horizon from temperature to temper- 
ature Ti ~ 1 GeV when the axion acquires mass. The decay of axion strings into axions has 



been discussed extensively in the literature and is reviewed in subsection [IIIB . 

At temperature Ti each string becomes the boundary of a domain wall. Walls not 
bounded by string are also present. We want to obtain statistical properties of this network 
of walls and strings under the assumption that the axion field is randomly oriented from one 
horizon to the next just before the axion mass turns on at time ti. In particular: What are 
the average densities of walls and of strings? Are there walls of infinite extent which are not 
cut up by any string? What is the energy fraction in walls not cut up by string? 

To address such questions, a cross-section of a finite but statistically significant volume 
of the universe near time ti was simulated in the following manner |Ty] . A 2-dim. hexagonal 



grid was constructed. Each hexagon represents a causally disconnected region. The circle 
{(0) = faC*" : — TT < a = ^ < +7r} is divided into three equal parts as shown in Fig. 2a. 
a = is at the middle of part 1. a = vr is at the boundary between parts 2 and 3. Each 
hexagon is then randomly assigned a number 1, 2 or 3. A hexagon to which #1 is assigned is 
thought of as a horizon volume in which the average value of a is in part 1 of the circle, and 
so on. Thus there is a domain wall at each boundary between two hexagons one of which 
has been assigned #2 and the other ^^3. Likewise there is an downward-going string at each 
vertex surrounded by hexagons to which #1, #2 and #3 have been assigned in clockwise 
order, and a upward-going string at each vertex surrounded by hexagons to which 7^1, 7^2 
and 7^3 have been assigned in counter-clockwise order. Fig. 2b shows some examples. Fig. 3 
shows a larger realization with the underlying grid removed. 
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Fig. 3 shows that there are no infinite domain walls which are not cut up by any string. 
The reason for this is simple. An extended domain wall has some probability to be cut up 
by string in each successive horizon it traverses. The probability that no string is encoun- 
tered after traveling a distance / along the wall decreases exponentially with /. Fig. 3 also 
shows that finite domain walls which are not cut up by string, and which therefore close 
onto themselves like a sphere or a do nut, are very rare. Indeed, the figure shows only one 
closed circle and a circle need not be a section of a sphere; it could be a section of a tube. 
The average density of walls at time ti predicted by the random process described here is 
approximately 2/3 per horizon volume. This is consistent with Fig. 3. 



III. THE COLD AXION POPULATIONS. 

In this section, we discuss the various contributions to the cosmological energy density 
in the form of cold axions. In addition to cold axions, there is a thermal axion population 
whose properties are derived by the usual application of statistical mechanics to the early, 
high temperature universe []T3[. The contribution of thermal axions to the cosmological 
energy density is subdominant if rua ^ 10^^ eV. We concern ourselves only with cold axions 
here. 

We distinguish the following sources of cold axions: 

1. vacuum realignment 

a. "zero momentum" mode 

b. higher momentum modes 

2. axion string decay 

3. axion domain wall decay 

The basis for distinguishing the contribution from vacuum realignment from the other two 
is that it is present for any quasi-Nambu-Goldstone field regardless of the topological struc- 
tures associated with that field. Contributions la and lb are distinguished by whether the 
misalignment of the axion field from the CP conserving direction is in modes of wavelength 
larger (la) or shorter (lb) than the horizon size ti at the moment the axion mass turns on. 
Contribution 2 is from the decay of axion strings before the QCD phase transition. Contri- 
bution 3 is from the decay of axion domain walls bounded by strings after the QCD phase 
transition. In the past literature on the axion cosmological energy density, contributions la 
and 2 have been discussed in detail whereas contributions lb and 3 have received much less 
attention. This section presents a systematic analysis of all contributions. 



A. Vacuum realignment 



Let us track the axion field from the PQ phase transition, when Upq{1) gets sponta- 
neously broken and the axion appears as a Nambu-Goldstone boson, to the QCD phase 
transition when the axion acquires mass. For clarity of presentation we consider in this 
subsection a large region which happens to be free of axion strings. Although such a region 
is rare in practice it may exist in principle. In it, the contribution to the cosmological axion 
energy density from string and domain wall decay vanishes but that from vacuum realign- 
ment does not. In more typical regions where strings are present, the contributions from 
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string and domain wall decay should simply be added to that from vacuum realignment. 
In the radiation dominated era under consideration, the space-time metric is given by: 

where t is the age of the universe, x are co-moving coordinates and R(t) ~ is the scale 
factor. The axion field a{x) satisfies: 

{d^D^ + ml{t)) a{x) = (d^ + 3^dt - -^V^ + ml{t)j a{x) = (3.2) 

where ma{t) = ma{T{t)) is the time-dependent axion mass. We will see that the axion mass 
is unimportant as long as ma{t) ^ 1/t, a condition which is satisfied throughout except at 
the very end when t approaches ti. The solution of Eq. (p.2|) is a linear superposition of 
eigenmodes with definite co-moving wavevector k: 

a{x,t) = J d^k a{k,t) e*^-^^ (3.3) 

where the a{k, t) satisfy: 

{d't+j^d, + ^^+mi{t)y{lt)=Q. (3.4) 

Eqs. ( ^.Ij ) and ( |3.3| ) show that the wavelength X{t) = ^R{t) of each mode is stretched by 
the Hubble expansion. There are two qualitatively different regimes in the evolution of each 
mode depending on whether its wavelength is outside (A(t) > t) or inside (A(t) < t) the 
horizon. 

For X{t) ^ t, only the first two terms in Eq. ( p. 4)) are important and the most general 
solution is: 

a{k, i) = ao{k) + a_i/2(A^)t"^/^ . (3.5) 

Thus, for wavelengths larger than the horizon, each mode goes to a constant; the axion field 
is so-called "frozen by causality" . 



For X{t) ^ t, the first three terms in Eq. (|3.4| ) are important. Let a{k, t) = R ^ {t)ip{k, t). 
Then 

(d''^+u\t))iP{k,t)={] (3.6) 

where 

-^W=™SW + ^ + ^-j|^. (3.7) 

Since uj -C c<j^, this regime is characterized by the adiabatic invariant 'ipQ{k^t)uj{t)^ where 
ipo{k,t) is the oscillation amplitude of ipik^t). Hence the most general solution is: 

a{K ^) = ^ cos (/ dt'uit')^ . (3.8) 
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The energy density and the number density behave respectively as g ~ -^3^ ~ -^^^ and 
n -. In -. Wttt, indicating that the number of axions in each mode is conserved. This 
is as expected because the expansion of the universe is adiabatic for X{t)t <^ 1. 

Let us define ^{uj,t) to be the number density, in physical and frequency space, of 
axions with wavelength A = — , for u > t~^. The axion number density in physical space is 
thus: 

n,(t)= / du;^{u;,t), (3.9) 
Ji-i aw 

whereas the axion energy density is: 

P dn 

Pa{t) = / duj-^{u,t)u . (3.10) 

Under the Hubble expansion axion energies redshift according to to' = uj^, and volume 

elements expand according to AV = AV (j^^ , whereas the number of axions is conserved 
mode by mode. Hence: 

dua, . ( R'Y dua R . . 

Moreover, the size of ^ for c<j ~ | is determined in order of magnitude by the fact that the 
axion field typically varies by fa from one horizon to the next. Thus: 



*"^^)f7V4(Vaf4|. (3.12) 



du^ ' ' „ . 1 du \t J \tj 2' ' 2 t2 

' t 

From Eqs. (|3.11|) and ( p.l2| ), and R ~ v^, we have: 

Eq. ( |3.13D holds until the moment the axion acquires mass during the QCD phase transition. 
The critical time is when ma{t) is of order t~^. Let us define ti. 

ma{h)h = 1 . (3.14) 

ma{T) was obtained from a calculation of the effects of QCD instantons at high temper- 
ature [§: 

. X q /lO^^GeVX /GeV\^ 

m,(T) ^ 4 ■ lO-'eV i j j (3.15) 

when T is near 1 GeV. The relation between T and t follows as usual from 

. / 1 \ ^ SttG 8nG , 
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where M is the effective number of thermal degrees of freedom. is changing near 1 GeV 
temperature from a value near 60, valid above the quark-hadron phase transition, to a value 
of order 30 below that transition. Using A/" — 60, one has 

t ~ 0.7 lO^O— — , 3.17 

sec Vsec/ \ fa 1 



which implies: 



The corresponding temperature is: 



■lO^^GeV^'/' 



Ti ^ IGeV y J j . (3.19) 

We will make the usual assumption 0] that the changes in the axion mass are adiabatic [i.e. 
m \t) ^it ^ ^a(^)] after ti and that the axion to entropy ratio is constant from ti till the 
present. Various ways in which this assumption may be violated are discussed in the papers 



of ref. |21|. We also assume that the axions do not convert to some other light axion-like 
particles as discussed in ref. . 

The above discussion neglects the non-linear terms associated with the self-couplings of 
the axion. In general, such non-linear terms may cause the higher momentum modes to 
become populated due to spinodal instabilities and parametric resonance. However, in the 



case of the axion field, such effects are negligible P,[23|l. 

We are now ready to discuss the vacuum realignment contributions to the cosmological 
axion energy density. 



1. Zero momentum mode 

This contribution is due to the fact that, at time ti, the axion field averaged over distances 
less than ti has in each horizon volume a random value between — vr/a and +TTfa, whereas 
the CP conserving, and minimum energy density, value is a = 0. The average energy density 
in this "zero momentum mode" at time ti is of order: 

pr'°(ii) - InilitM ■ (3.20) 

The corresponding average axion number density is: 

1 f2 

[tl) = 7-Tpa (^l) ~ TTT • l^-^lj 

ma{ti) 2ti 

Since ma{t) is assumed to change adiabatically after ti, the number of axions is conserved 
after that time. Hence: 




(3.22) 



The axions in this population are non-relativistic and contribute each to the energy. 
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2. Higher momentum modes 



This contribution is due to the fact that the axion field has wiggles about its average 
value inside each horizon volume. The axion number density and spectrum associated with 
these wiggles is given in Eq. (|3.13|) . Integrating over oj > t^^, we find: 

nr'\t) = nZ-\t,) (^)' ~ 1| (^)' (3.23) 

for the contribution from vacuum realignment associated with higher momentum modes. 
The bulk of these axions are non-relativistic after time ti and hence each axion contributes 
nia to the energy. Note that n^"'^'°(t) and n^"'^'^(t) have the same order of magnitude. 



B. String decay 

The contribution from axions which were produced in the decay of axion strings has been 
discussed extensively in the literature. We describe it here for the purpose of completeness 
and clarity. There is still disagreement about the size of this contribution. We will go rapidly 
over those points which are not controversial, but indicate explicitly those which are. 

Axion strings have energy per unit length: 

~ r dk^ = 7ivl\n{vaL) (3.24) 
JL-i k 

where L is an infra-red cutoff. This is the same formula as Eq. ( p.3|) but with the integral 
written in wavevector space. In the early universe, when the strings have random directions 
relative to one another, L is of order the distance between neighboring strings. The strings 
move relativistically and decay efficiently into axions. As a result there is approximately 
one long string per horizon at all times after the temperature has dropped below T^, but 
before the QCD phase transition. By definition, a long string traverses the horizon. By 
intersecting each other with reconnection, the long strings produce closed string loops of 
size smaller than the horizon which decay efficiently into axions. 
One long axion string per horizon implies the energy density: 

PsUt) ^ -^\n{tva) . (3.25) 

We are interested in the number density of axions produced in the decay of axion strings 
because, as we will soon see, these axions become non-relativistic soon after time ti and 
hence contribute each m^, to the energy. The equations governing the number density n^*^ (t) 
of axions produced in the decay of strings are: 

^ = -2Hp,tr - (3.26) 
dt ^ dt ^ ' 



dn 



str 



dt 



-3HnT' + 



1 dp. 



str—*a 



00 (t) dt 



(3.27) 
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where uj{t) is defined by: 



(3.28) 



In Eqs. (13:261^:281) , ^^^^^{t) is the rate at which energy density gets converted from strings 



into axions at time t, whereas is the spectrum of axions thus emitted. The term 

—2Hpstr = +Hpstr — SHpstr in Eq. (|3.26| ) takes account of the fact that the Hubble ex- 
pansion both stretches {+Hpstr) and dilutes {—3Hpstr) long strings. To obtain n'^a^it) from 
Eqs. ( |3.25[ - p:^BD , we must know what uj{t) characterizes the spectrum l^t^^" of axions 
radiated by cosmic axion strings at time t. 

The axions that are radiated at time t are emitted by strings which are bent over a 
distance scale of order t and are relaxing to lower energy configurations. The main source 
is closed loops of size t. Two views have been put forth concerning the motion and the 
radiation spectrum of such loops. One view is that such a loop oscillates many times, 



with period of order t, before it has released its excess energy and that the spectrum of 
radiated axions is concentrated near Let us call this scenario A. A second view PJT2 



is that the loop releases its excess energy very quickly and that the spectrum of radiated 
axions '^^f'^" ~ ^ with a high energy cutoff of order Va and a low energy cutoff of order 
Let us call this scenario B. In scenario A one has uj{t) ~ hence, for t < ti. 



nf'^{t)^'^Hvat) , (3.29) 



t 

2 



whereas in B one has ti;(t) ~ ^ \n(vat) and hence 



2 



nf '^(t) ~ ^ . (3.30) 



Comparing Eqs. (|3l3|) , ( ^291) and (|3o]), one sees that in A the contribution from string 



decay to the cosmological axion energy density dominates by a factor \n(vat) ~ 100 over the 
contribution from vacuum realignment, whereas in B the two contributions are of the same 
order of magnitude. 

We have carried out computer simulations ||12| of the motion and decay of axion strings 
with the purpose of deciding between the two possibilities. In scenario A, the wavevector 
spectrum ^ of the energy stored in the axion field changes from ^Ji^Hit ^ characteristic of the 
string [see Eq. ( ^.24| )], to a spectrum peaked near ~ after the string loop has decayed 



into axions. In scenario B, ^ does not qualitatively change. To distinguish between the 
two possibilities, we defined the quantity: 

N^^jdk--, (3.31) 

and computed it as a function of time while integrating the equations of motion which follow 
from Lagrangian (|2.1| ) with 7] = In scenario A, Nax increases as the spectrum softens, 
whereas in 5, Nax remains approximately constant. We found in our simulations that Nax 
decreases by about 20% during the decay of closed string loops. In the simulations, In(t'aL) 
is in the range 4 to 6. If scenario A were correct, A^^^. should have increased by the factor 
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In(faL) — 1 = 3 to 5 (depending upon the loop size L) during the decay of loops. So our 
simulations support scenario B. Note that Nax is the quantity that directly measures the 
string contribution to the axion cosmological energy density. An important caveat is that 
the simulations are done for values of \a{vaL) much smaller than those, of order 60, which 
characterize axion strings in the early universe. 

We will assume henceforth that axion production by cosmic string decay is consistent 
with possibility B. Hence: 

nTit)-f{^)^ ■ (3-32) 



ti\R 



It is of the same order of magnitude as the contributions [Eqs. (|3.22| ) and (|3.23|) ] from 
vacuum realignment. 



C. Domain wall decay 

During the QCD phase transition, when the axion acquires mass, the axion strings 
become the boundaries of axion domain walls. Fig. 3 shows a cross-section of a piece of the 
universe at that time. In 3 dimensions, the domain walls bounded by string may be like 
pancakes or they may be long surfaces with holes. Which type of configuration dominates 
matters little to the arguments given below. It does matter that the universe is relatively 
free of domain walls which are not bounded by strings, i.e. domain walls which close onto 
themselves like spheres, donuts and so on. We concluded in section |I| that such domain 
walls are in fact very rare. 

Near 1 GeV temperature, the time- dependent axion mass is given approximately by 



Eq. ( |3.17|) . The wall energy per unit surface is 

a{t) 8/X(t). (3.33) 

Since the chiral symmetry breaking phase transition has presumably not taken place yet, 
we do not use the zero-temperature result Eq. (|2.6|). As illustrated in Fig. 1, the string at 
the boundary of a wall is embedded into the wall. Hence its infra-red cutoff L, in the sense 
of Eqs. ( p.3|) and ( p.24| ), is of order and its energy per unit length is therefore 



/i~7r/2ln(/,/mJ . (3.34) 

The surface energy of a typical (size ~ ti) piece of wall bounded by string is cr(t)t1 
whereas the energy in the boundary is ~ fiti. There is a critical time ^2 when the ratio 

E^ Tcln{fa/ma) 
is of order one. Using Eqs. ( |3.17D and ( ^.18| ), we (crudely) estimate : 
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(3.37) 



After t2 the dynamics of the walls bounded by string is dominated by the energy in the walls 
whereas before ^2 it is dominated by the energy in the strings. A string attached to a wall 
is pulled by the wall's tension. For a straight string and flat wall, the acceleration is: 

aJt) = ~ , ,"V , ~ ~ 3.38 

/i 7rln(/„/m„) 23 ti m,(t2) ^ ^ 

Therefore, after t2, each string typically accelerates to relativistic speeds, in the direction of 
the wall to which is it attached, in less than a Hubble time. We argue below that the string 
will then unzip the wall, and that the energy stored in the wall gets released in the form 
of barely relativistic axions. However, before doing so, let us describe the competing decay 
mechanisms which have been discussed. 

One mechanism is emission of gravitational waves . A piece of wall bounded by string 
of size i, assuming for the moment that it is long lived and that its motion is not damped by 
friction against the surrounding plasma, oscillates with frequency ~ and hence emits 
gravitational waves. The power P may be estimated using the quadrupole formula: 

P ~ ~ Giatfuj'' ~ Ga^e^ . (3.39) 

Indeed, emission of gravitational waves by walls bounded by string of size ti is dominated by 
the wall contribution, rather than the string contribution, after time t2- Eq- (3.39) implies 
the lifetime: 

Vav~(Ga)-i:^3.10='sec(^^^) , (3.40) 

independently of size. Note that a piece of wall bounded by string may self-intersect with 
reconnection producing two walls bounded by string which may again intersect with recon- 
nection, and so on. This fragmentation does not affect the lifetime ( |3.4CI| ) although it shifts 
the spectrum of gravitational waves to higher frequencies. 

Another energy dissipation mechanism is drag caused by the reflection of particles in the 
primordial plasma. The reflection coefficients between various particles and axion domain 
walls are discussed in ref. It was found there that axions which are not highly relativistic 
have reflection coefficient of order one and that the dominant friction experienced by the 
oscillating walls is on the fluid of cold axions that were produced by vacuum realignment 
and string decay. Consider a bent wall of curvature radius of order ti. The tension in the 
wall provides a driving force per unit surface of order a/ti whereas reflection of cold axions 
produces a damping force per unit surface of order 2mal3'^na where (3 is the velocity of the 
wall. Adding the contributions from vacuum realignment [Eqs. ( p.21j ) and (|3.23|) ] and from 

string decay [Eq. (|3.32| )] we have Uait) ~ ^ . This neglects the contribution n^"'(t) 

of axions produced in the decay of the walls themselves but we will see below that n'^^it) IS 
of the same order of magnitude as n^"'^(t) and n'^J^^{t). Comparing the driving force with the 
damping force, one concludes that (3 is of order one immediately after ti and hence that the 
rate of energy dissipation by friction is P ~ — '^'■^^ ^ ~ mana{t)l'^ which implies: 
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dt ti\t 

Thus energy dissipation by friction is important just after the domain walls appear but this 
mechanism soon turns off as the universe becomes dilute. 

If emission of gravitational radiation and friction on the cold axion fluid were the only 
important dissipation mechanisms, the conclusion would be that a typical wall, whether or 
not bounded by string, loses a large fraction of its energy by friction immediately after ti 
and then lives till a time of order Tgrav when it decays into gravitons. 

However we find it far more likely that the walls bounded by string decay into axions, long 
before Tgrav If emission of gravitational waves were the dominant dissipation mechanism, 
a piece of wall bounded by string of size i, assuming for the moment that it does not 
reduce its size by self-intersections with reconnection, would oscillate on the order of ^^j^ ~ 

10^° (^Y^ ( io"^5°cv ) ^ times before decaying away. This seems an implausibly large number 
considering that nothing forbids the wall from decaying into axions instead. It is especially 
implausible when one considers that at each oscillation a wall bounded by string may self- 
intersect with reconnection producing two smaller walls bounded by string. If the two pieces 
produced by an intersection with reconnection are typically of roughly equal size it takes 
only log2(mati) ~ 10 iterations of this process for the size of the pieces to go down from ti 
to the thickness of the walls. Pieces of size surely decay into axions because their 
dynamics is dominated by the energy of the strings and strings are known to decay into 
axions efficiently. Even if the two pieces produced by a self-intersection with reconnection 
are typically of very different size, it takes at most tiiria iterations for the pieces to go 
down in size to m~^. Hence for emission of gravitational waves to be the dominant decay 
mechanism, the probability p of self-intersection with reconnection per oscillation can not 

be larger than ^ 10^^ ^ev \ / ^ reasonably be thought to be of order one. It 

'^grav \ '^a / 

seems implausible that p could be as small as 10~^. 

In section |IV| we describe our computer simulations of the motion and decay of walls 
bounded by string. In the simulations the walls decay immediately, i.e. in a time of or- 
der their size divided by the speed of light. The average energy of the radiated axions is 
(ua) ~ in the simulations. However, as was already mentioned in the Introduction, 
the simulations are done for ln{va/ma) — In(lOO) ^ 4.6, whereas in axion models of interest 
ln{va/ma) — 60. One reason this may affect the results is the following. When massless 
the axion field is equivalent, by a duality transformation, to an anti-symmetric two index 
gauge field A^^^, [^. The axion string is a source for this field. If the string accelerates 
with acceleration it emits radiation of frequency of order a^. Since ~ ^ in axion 
models of interest [see Eq. ( ^.381 )], the emission of axion radiation is inhibited by the fact 



that as is considerably less than rria, the lowest possible frequency for axion radiation. In 
the simulations as — 0.6 ma. 

At any rate, for the reasons given earlier, it is most plausible that the walls bounded by 
string decay into barely relativistic axions long before Tgrav Let be the time when the 
decay effectively takes place and let 7 = ^^^^^ be the average Lorentz 7 factor then of the 
axions produced. In section II, we estimated the density of walls at time ti to be of order 0.7 
per horizon volume. Hence we estimate that between ti and the average energy density 
in walls is 
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Pd.w.W~(0.7)(9)m.(t)^(:|)'. (3.42) 

We used Eq. ( ^.61 ) and assumed that the energy in walls simply scales as niait). After time 
ts, the number density of axions produced in the decay of walls bounded by strings is of 
order 

^d.w.^^) ^ (3.43) 



{iUa) \RJ -fh\R 

Note that the dependence on drops out of our estimate of n"^-^-. If we use the value 7 = 7 
observed in our computer simulations, we have 

ni-it) ~ f (f )' • (3-44) 



In that case the contributions from vacuum realignment [Eq. ( p.22|) and Eq. ( p. 231) ] and 



domain wall decay are of the same order of magnitude. However, we will find in section 
IV that 7 rises approximately linearly with \n{\/Xva/ma) over the range of \n{y/Xva/ma) 
investigated in our numerical simulations. If this behaviour is extrapolated all the way to 
\n{\/Xva/ma) — 60, which is the value in axion models of interest, then 7 ~ 60. In that case 
the contribution from wall decay is subdominant relative to that from vacuum realignment. 



D. The cold axion cosmological energy density 

Adding the three contributions, and assuming scenario B for the string contribution, we 
estimate the present cosmological energy density in cold axions to be: 

Paito)r^3^(^)'ma. (3.45) 



ti \Ro 

Following 1^, we may determine the ratio of scale factors Ri/Rq by assuming the conser- 
vation of entropy from time ti till the present. The number of effective thermal degrees of 
freedom at time ti is approximately A/i — 61.75. We do not include axions in this number 
because they are decoupled by then. Let ^4 be a time (say T4 = 4 MeV) after the pions 
and muons annihilated but before neutrinos decoupled. The number of effective thermal 
degrees of freedom at time is A/4 = 10.75 with electrons, photons and three species of 
neutrinos contributing. Conservation of entropy implies N'lTfRf = J^f^T^R^. The neutrinos 
decouple before e~^e~ annihilation. Therefore, as is well known, the present temperature 
Tyo = 2.735 ''K of the cosmic microwave background is related to T4 by: ^^1^1 = 2 Rl 
Putting everything together we have 

p \ 3 ,rTi ^ 3 



-^j ^ 0.063 (^^) . (3.46) 



Combining Eqs. (|3A8D , and ( CT ), 



.-29 / fa 



\ 7/6 
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Dividing by the critical density pc = we find : 

^ fa y^'(o.7-' 

^2GeV J \ h 

where h is defined as usual hj Hq = h lOOkm/s-Mpc. 



E. Pop. I and Pop. II axions 

On the basis of the discussion in section IIIA-D we distinguish two kinds of cold axions: 

I) axions which were produced by vacuum realignment or string decay and which were 
not hit by moving domain walls. They have typical momentum (piUi)) ~ 7- at time ti 
because they are associated with axion field configurations which are inhomogeneous on the 
horizon scale at that time. Their velocity dispersion is of order: 




Let us call these axions population I. 

II) axions which were produced in the decay of domain walls. They have typical mo- 
mentum (^//(ts)) ~ 7771^(^3) at time when the walls effectively decay. Their velocity 
dispersion is of order: 

Mt) ~ . 10- ( ,^^] ('-^\'° ^ . (3.50) 

rUa R \ TUa RlJ \ rUa J R 

The factor q = 1 ^^^^^ "il P^i-rameterizes our ignorance of the wall decay process. We expect 
q to be of order one but with very large uncertainties. Fortunately there is a lower bound on 
q which follows from the fact that the time ^3 when the walls effectively decay must be after 
t2 when the energy density in walls starts to exceed the energy density in strings. Using 
Eq. ( g:36D , we find: 



^ _ ^majts) Rs ^ im,{h) R " /in-5„..x2/3 



1 130 V m„ y ^ ' 



rUa Ri nia R 

Since our computer simulations suggest 7 is in the range 7 to 60, the axions of the second 
population (pop. II) have much larger velocity dispersion than the pop. I axions. 

Note that there are axions which were produced by vacuum realignment or string decay 
but were hit by relativistically moving walls at some time between ti and ^3. These axions 
are relativistic just after getting hit and therefore are part of pop. II rather than pop. I. 

The fact that there are two populations of cold axions, with widely differing velocity 
dispersion, has interesting implications for the formation and evolution of axion miniclusters. 
This is discussed in section 0. 
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Here we speculate that the primordial velocity dispersion of pop. II axions may some 
day be measured in a cavity- type axion dark matter detector [Q. If a signal is found in 
such a detector, the energy spectrum of dark matter axions will be measured with great 



resolution. It has been pointed out that there are peaks in the spectrum ||T^ because late 
infall produces distinct flows, each with a characteristic local velocity vector. These peaks 
are broadened by the primordial velocity dispersion, given in Eqs. ( |3.49|) and (|3.5CI| ) for pop. I 
and pop. II axions respectively. These equations give the velocity dispersion in intergalactic 
space. When the axions fall onto the galaxy their density increases by a factor of order 10^ 
and hence, by Liouville's theorem, their velocity dispersion increases by a factor of order 10. 
[Note that this increase is not isotropic in velocity space. Typically the velocity dispersion 
is reduced in the direction longitudinal to the flow in the rest frame of the galaxy whereas 
it is increased in the two transverse directions.] The energy dispersion measured on Earth 
is /\E = rriaPAf] where (3 ~ 10~^ is the flow velocity in the rest frame of the Earth. Hence 

-19 f 10-5eV A^/Q_ fr.r.^r.^ To^^^v^oov.^ A?? 10"^^- / lO-'-^eV \ 



we find AEj ~ 3 ■ 10 Ao ^ey \ pop. I axions and AEjj 



ma 



for pop. II. The minimum time required to measure AE is (AE)"^. This assumes ideal 
measurements and also that all sources of jitter in the signal not due to primordial velocity 
dispersion can be understood. There is little hope of measuring the primordial velocity 

10-5eV\l/6 



dispersion of pop. I axions since {AEj 

■lO-5eV\5/6 



10 years (i^;^) • However {AEn] 



-1 



day 



q ^ ; and hence it is conceivable that the primordial velocity dispersion of pop. II 



axions will be measured. 



IV. COMPUTER SIMULATIONS 



We have carried out an extensive program of 2D numerical simulations of domain walls 
bounded by strings. The Lagrangian (2.1) in finite difference form is 



0i(^ + - M^^t) 



(4.1) 



where n labels the sites. In these units, Va = 1, the wall thickness is 1/ma = 1/^^ the 
core size is 5 = 1/a/A. The lattice constant is the unit of length. In the continuum limit, 
the dynamics depends upon a single critical parameter, namely rUaS = rria/ \fX. 

Large two-dimensional grids (~ 4000 x 4000) were initialized with a straight or curved 
domain wall, at rest or with angular momentum. The initial domain wall was obtained 
by overrelaxation, starting with the Sine-Gordon ansatz 0i + i02 = exp(42 tan~^exp(ma|/)) 
inside a strip of length the distance D between the string and anti-string, and the true 
vacuum (0i = 1, 02 = 0) outside. The string and anti-string cores were approximated by 
(pi + i(f)2 = — tanh(.583r/5) exp(=Fi6'), where r and 6 are polar coordinates about the core 
center, and held fixed during relaxation. Stable domain walls were obtained for l/(ma5) ^ 
3. A first-order in time and second-order in space algorithm was used for the dynamical 
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evolution, with time step dt = 0.2. The boundary conditions were periodic throughout and 
the total energy was conserved to better than 1%. If the angular momentum was nonzero, 
the initial (p was obtained as the difference of two relaxed wall configurations a small time 
step apart, divided by that time step. 

The evolution of the domain wall was studied for various values of VX/rUa, the initial 
wall length D and the initial velocity v of the strings in the direction transverse to the wall, 
the string and anti-string going in opposite directions. The strings attached to the wall are 
rapidly accelerated by the wall tension, the potential and gradient energies of the wall being 
converted into kinetic energy of the strings. Fig. 4 shows the longitudinal velocity of the 
core as a function of time for the case — 400, y/X/nia — 10 and v — 0. An important 
feature is the Lorentz contraction of the core. For reduced core sizes S/'js ^ 5, where 7^ 
is the Lorentz factor associated with the speed of the string core, there is 'scraping' of the 
core on the lattice accompanied by emission of spurious high frequency radiation. This 
artificial friction eventually balances the wall tension and leads to a terminal velocity. In 
our simulations we always ensured being in the continuum limit. 

For a domain wall without rotation {v = 0), the string cores meet head-on and go through 
each other. Several oscillations of decreasing magnitude generally occur before annihilation. 
For 73 ~ 1 the string and anti-string coalesce and annihilate one another. For 7^ ^ 2, the 
strings go through each other and regenerate a new wall of reduced length. The relative 
oscillation amplitude decreases with decreasing coUision velocity. Fig. 5 shows the core 
position as a function of time for the case m~ ^ = 1000,v^/m„ = 14.3, D = 2896. Fig. 6a 
shows the gradient, kinetic, string potential, wall potential and total energies as a function 
of time. By 'string (wall) potential' energy, we mean the energy associated with the third 
(fourth) term in Eq. (4.1). The total energy is accurately conserved in spite of the violent 
goings-on when the string and anti-string meet. 

We also investigated the more generic case of a domain wall with angular momentum. 
The strings are similarly accelerated by the wall but string and anti-string cores rotate 
around each other. Fig. 6b shows the gradient, kinetic, string potential, wall potential and 
total energies as a function of time for the case m~ 1 = 500, y/X/nia = 28.3, D = 2096, and 
V — 0.25. The field is displayed in Fig. 7 at various time steps for the case — 100, A = 
0.01, D = 524 and v = 0.6. No oscillation is observed. All energy is converted into axion 
radiation during a single collapse. Oscillations occur only when the angular momentum is 
too small for the string and anti-string cores to miss each other. 

We performed spectrum analysis of the energy stored in the field as a function of time 
using standard Fourier techniques. The two-dimensional Fourier transform is defined by 



fip) = ^==Eexp 



\ Ly ^ 



fin) (4.2) 



for = 1, • • • , Lj where j — x,y. The dispersion law is: 



COr, = 



p - \ 2 2 - cos — cos — — ^ + 171^. (4.3) 



Fig. 8 shows the power spectrum dE /dloguj of the field at various times during the decay 
of a rotating domain wall for the case — 500, A = 0.0032, D — 2096, and v — 0.25. 
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Initially, the spectrum is dominated by small wavevectors, k ~ ma. Such a spectrum is 
characteristic of the domain wall. As the wall accelerates the string, the spectrum hardens 
until it becomes roughly 1/k with a long wavelength cutoff of order the wall thickness 1 / 
and a short wavelength cutoff of order the reduced core size S/'-ys- Such a spectrum is 
characteristic of the moving string. In the case of the figure, the string and anti-string 
cores annihilate at about the time of the fourth frame (Fig. 8d). The spectrum remains 
approximately 1/k during the annihilation process and thereafter. 

As discussed in section IIIC, 7 = {u)a)/ma, the average energy of the radiated axions 
in units of the axion mass, is the quantity which determines the wall decay contribution to 
the axion cosmological energy density [see Eq. (3.43)]. Fig. 9 shows the time evolution of 
{uj)/ma for = 500, = 0.25, D = 2096 and various values of A. By definition, 

H = E^p/E- (4.4) 

p p y 

where is the gradient and kinetic energy stored in mode p of the field 0. After the domain 
wall has decayed into axions, {u) = {ua). Fig. 10 shows the time evolution of {u)/ma for 
= 500, A = 0.0016, D = 2096 and various values of v. For small v and/or small A, i.e 
when the angular momentum is low and the core size is big, the string cores meet head on 
and the strings have one or more oscillations. In that case, {uja)/ma — 4. This is consistent 



with the value, {uja)/ma — 3, found by Nagasawa and Kawasaki |jT^. The low angular 
momentum regime, when the strings oscillate, is characterized by low energy of the radiated 
axions. We believe this regime to be less relevant to wall decay in the early universe because 
it seem unlikely that the angular momentum of a wall at the QCD epoch could be small 
enough for the strings to oscillate then, especially when one considers that the actual wall 
configurations are 3 dimensional rather than 2D. 

In the more generic case when no oscillations occur, {ua) / ma — 7. Moreover, we find 
that {uja)/ma depends upon the critical parameter \/X/ma, increasing approximately as the 
logarithm of that quantity; see Fig. 9. This is consistent with the time evolution of the energy 
spectrum, described in Fig. 8. For a domain wall, (u) ~ whereas for a moving string 
(u) ~ ma \n{^/Xva'ys/ma). Since we find the decay to proceed in two steps: 1) the wall energy 
is converted into string kinetic energy, and 2) the strings annihilate without qualitative 
change in the spectrum, the average energy of radiated axions is (uja) ~ ln(^/A^'a/'"^■a)■ 
Assuming this is a correct description of the decay process for ^/Xva/ma ~ 10^^, then 
{u!a)/ma ~ 60 in axion models of interest. 

Present technology is inadequate for 3D simulations with sufficient resolution for the 
purposes we are interested in, such as obtaining reliable estimates of the factor 7. However, 



lower resolution 3D and 2D simulations of string-wall networks were carried out in ref. ^6 
They showed qualitatively similar behaviour in 2 and 3 dimensions. 



V. AXION MINICLUSTERS 

We saw in section III that the axion fluid is inhomogeneous with ^ = 0{1) at the 
time of the QCD phase transition. As will be discussed below, the streaming length of the 
pop. I axions (those from vacuum realignment and string decay) is too short for all these 
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inhomogeneities to get erased by free streaming before the time teq of equality between matter 
and radiation, teg is when density perturbations start to grow by gravitational instability. 
At time tgq, the ^ = 0{1) inhomogeneities in the axion fluid promptly form gravitationally 



bound objects, called axion miniclusters |ll6| , p!7| . Of course, axion miniclusters occur only if 



there is no inflation or if inflation occurs with reheat temperature above the phase transition 
where Upq{1) is spontaneously broken. 

Axion miniclusters were discussed in the seminal papers of Hogan and Rees []T6[ , and of 



Kolb and Tkachev [|1^. However, in these pioneering studies the role of domain walls, and 
hence of axions from domain wall decay, was neglected. Also the estimates of the minicluster 
mass have varied considerably from one another. Hence we think it appropriate to give our 
own analysis of this interesting phenomenon. 

As described in section HI, there are two populations of cold axions, pop. I and pop. II, 
with velocity dispersions given by Eqs.( ^.41^ ) and ( p. 501 ) respectively. Both populations are 
inhomogeneous with ^ = 0{1) at the time of the QCD phase transition. The free streaming 
length from time ti to teq is: 

..OU^ IT^.^fteq 



= R{Ug) J^^ dt^^ ^ /3(ti)^tite,ln (^^j . (5.1) 

For the sake of definiteness, we assume henceforth Qq = 1 and Hq = 60 km/s.Mpc, in which 
case Teq = 2 eV. The free streaming length should be compared with the size 

i^e -h^^ \ft^q ^ 0.7 ■ lO^^cm (5.2) 
of axion inhomogeneities at time t^q. Using Eq. (|3.49|) we find for pop. I: 



-^,4^)^1.6. 10- fi2:!!Xy". (5.3) 

txTTia \ti J \ rUa J 



Hence, in the axion mass range of interest, pop. I axions do not homogenize. At teg many 
pop. I axions condense into miniclusters since they have ^ = 0{1) upon arrival. The typical 
size of axion miniclusters is imc and their typical mass is: 



M^, ~ ^ Pa{teq)f J ~ 0.7 ■ IO'^^Mq (^^^^^ • (5-4) 

We assumed that half the cold axions are pop. I and that half of those form miniclusters, 
hence the factor 1/4. We used Eq. (3.48) to estimate Paiteq)- 
Using Eq. (|3.50| ), we find for pop. II: 



glnl ^) ~40g. (5.5) 

Using Eq. (3.51) and the range 7 ~ 7 to 60, suggested by our numerical simulations, we 
conclude that pop. II axions do homogenize and hence that the axion energy density has a 
smooth component at t^q. 
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At time t, pop. II axions get gravitationally bound to the nearest minicluster if it is less 
than distance d ~ GMmd away. Since most are at a distance d{t) ~ imc^fj, the bulk 
of pop. II axions become gravitationally bound to miniclusters when the temperature has 
dropped to: 



'II 



4-10-=^%, (5.6) 



assuming ~ 1. Thus we arrive at a qualitative picture of miniclusters before galaxy 
formation as having an inner core of pop. I axions and a fluffy outer envelope of pop. II 
axions. The inner core has size f^c and density 

Pmc ~ 10-i«^ . (5.7) 
cm'^ 

If Tji > To (i.e. g ^ 6), the outer envelope has roughly the same mass as the core but size 

2 lO^^cm (5.8) 
and density 

Pmc- 10-^^^ . (5.9) 

If Tjj < To, then the envelope is correspondingly less massive and less dense. When a 
minicluster falls onto a galaxy, tidal forces of the galaxy are apt to destroy it. We find below 
that the outer envelopes almost surely get pulled off whereas the inner cores barely survive. 

When a minicluster passes by an object of mass M with impact parameter b and velocity 
V, the internal energy per unit mass AE given to the minicluster by the tidal gravitational 



forces from that object is of order 116 



whereas the binding energy per unit mass of the minicluster E G Pmc^mc- the minicluster 
travels a length I = (3t through a region where objects of mass M have density n, the relative 
increase in internal energy is: 

(5.11) 

^ Pmc 



where pm = Mn. Eq. ( |5.11| ) follows from the fact that AE is dominated by the closest 
encounter and the latter has impact parameter hmin'- '^^'inin^^ = 0(1)- Note that ^ is 
independent of M. The effect upon a minicluster of falling once with velocity [3 ~ 10~^ 
through the inner halo (r < 10 kpc) of our galaxy where Pm ~ 10~^^gr/cm^ is: 

P^' in-6 f +1. • 
— - ~ ~ 10 tor the mner core 

E Pmc 

~ 10 for the outer envelope (5.12) 
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Since we expect q to be of order one, it seems very likely that the envelope gets pulled off 
after one or more crossings of the inner parts of the galaxy. This result is reassuring for 
direct searches of dark matter axions on Earth since it implies that a smooth component of 
dark matter axions with density of order phaio permeates us whether or not there is inflation 
after the Peccei-Quinn phase transition. If there were no pop. II axions one might fear that 
nearly all axions are in axion miniclusters. 

This result can also be obtained by considering the tidal radius of a mini-cluster at our 
location in the galaxy: 

Rtidai ~ {inr^f'^D ~ 8 • lO^^cm 5.13 

^Mgal \ rUa J 

where D ~ 8.5 kpc is our distance from the galactic center and Mgai — IO^^Mq is the mass 
of the galaxy interior to us. Rudai is larger than the inner core size imc but smaller than the 
outer envelope i'^^. 

A minicluster inner core which has spent most of its life in the central part of our galaxy 
only barely survived since ^ ~ 10~^ in that case. If most minicluster inner cores have 
survived more or less intact, the direct encounter of a minicluster with Earth would still be 
quite rare, happening only every 10'^ years or so. The encounter would last for about 3 days 
during which the local axion density would increase by a factor of order 10^. 



VI. THE THIRD SOLUTION TO THE AXION DOMAIN WALL PROBLEM 

In this section we consider axion models in which the number of degenerate vacua at 
the bottom of the Mexican hat potential is larger than one. N is an integer given by 

Tr(gpQg^g^) = ns'^". (e.i) 

Here Qpq is the Peccei-Quinn charge, (a = 1...8) are the color charges and Tr represents 
a sum over all the left-handed Weyl spinors in the model. It is easy to construct models for 
an arbitrary value of N. 

There are N distinct degenerate equally spaced vacua at the bottom of the Mexican hat 
potential |^. The vacua are related by an exact spontaneously broken Z{N) symmetry. A 
domain wall is a transition in space from any one vacuum to a distinct neighboring vacuum. 
A string is the boundary of A^ walls. If A^ > 1 and there is no inflation after the Peccei- 
Quinn phase transition, axion models suffer from a domain wall problem because the energy 
density in walls ends up exceeding the energy density in matter and radiation. A domain 
wall dominated universe is very much unlike our own. The walls are gravitationally repulsive 



25| . The universe is divided into rapidly expanding bubbles with walls at the boundaries 
and islands of matter in the middles. When averaged on scales larger than the bubble size, 
the universe expands as i? ~ At least in the case of axion walls, the islands are far 



too small for one island to be identified with our visible universe. 

The axion domain wall problem can be solved 1) by having inflation with reheating 
temperature less than the PQ phase transition temperature, or 2) by having A^ = 1 as 
discussed in the other sections of this paper, or 3) by having a small explicit breaking of the 
Z{N) symmetry [Q. The third solution is discussed in this section. The symmetry breaking 
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must lift completely the degeneracy of the vacuum and be large enough that the unique true 
vacuum takes over before the walls dominate the energy density. On the other hand, it must 
be small enough that the PQ mechanism still works. The implications of biased discrete 
symmetry breaking have been discussed in a general manner in ref. ||2^ . 



We may, for our purposes, take the effective action density for the axion field to be: 



1 „ mlv'^ 



Ca = -daad^a + 



a" a 



2 ^ iV2 



cos 



(Na 

V Va 



(6.2) 



The axion field a{x) is cyclic with period 2nva- Thus Eq. ( |672D implies N degenerate vacua. 
The axion decay constant is /„ = The domain wall is a transition in space from a 
vacuum, say a = 0, to a neighbor vacuum, say a = The energy per unit surface is 

cr = Sniafa- The actual axion domain wall is only crudely described by Eq. (|6.2| ) because it 
involves not only a rotation of the axion field but also of the quark-antiquark condensates. 
However Eq. (6.2) is adequate for the order of magnitude estimates we are interested in here. 

Between the PQ and QCD phase transitions, axion strings are present and evolve as in 
the = 1 case. At time ti each string becomes the boundary of N walls. For a brief period 
the walls lose energy by friction against the primordial plasma but this dissipation process 
is soon turned off as the plasma becomes dilute. 

Causality implies that there is at least one domain wall per horizon on average since each 
horizon volume picks a vacuum independently from its neighbor. If there is of order one 
wall per horizon, the energy density in walls p ~ ^ = j. In that case, the energy in walls 
per co-moving volume {E ~ pR^ ~ ^/t) increases with time. This suggests that there is not 
more than approximately one wall per horizon because energy in walls has low entropy and 
hence increases only reluctantly. 

With one wall per horizon, the fraction of the energy density in walls is 

fi....(t)~^^^~2.10-«(^:^) -. (6.3) 
3 \ rua J h 

If this fraction reaches one, the universe becomes domain wall dominated and disaster strikes 
as described above. To get rid of the walls, we add to the model a tiny Z{N) breaking term 
which hfts the vacuum degeneracy completely, e.g.: 

SV = + h.c.) (6.4) 

where (p is the PQ field with the "Mexican hat" potential, at the bottom of which '^{x) = 
Va e*"*^^-'/^" . On the RHS of Eq.(6.2) an extra term appears: 

6Ca = 2va^ COS ( - - s) . (6.5) 

Now the unique true vacuum is the one for which 15 — ^ I is smallest. Its energy density is 
lowered by an amount of order relative to the other quasi- vacua. As a result, the walls at 
the boundary of a region in the true vacuum are subjected to an outward pressure of order 
^Va- Since the walls are typically a distance t apart, the volume energy ^ Va associated 
with the lifting of the vacuum degeneracy grows more rapidly than the energy at'^ in the 
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walls. At a time r ~ the pressure favoring the true vacuum starts to dominate the 
wall dynamics and the true vacuum takes over, i.e. the walls disappear. We require this 
to happen before the walls dominate the energy density. Using Eqs. (|3.18| ) and ( |6.3| ), we 
obtain: 

r ~ < 10^ sec (-^] . (6.6) 



^Va VlO-^eV. 
On the other hand ^ is bounded from above by the requirement that SV does not upset 
the PQ mechanism. 6V shifts the minimum of the effective potential for the axion field, 



inducing a 6^ ~ — li^- The requirement that 9 < 10 implies: 



r>.ec('^]. (6.7) 

Eqs. ( |6.6| ) and ( |6.7|) indicate that there is very little room in parameter space for this third 
solution to the axion domain wall problem, at least if rria ~ 10~^ eV. But the third solution 
is not completely ruled out. 

Soon after time r, a cross-section of a small fraction of the universe looks like Fig. 3 
except that every string - anti-string pair is connected by N walls instead of just one wall. 
This likely changes the mechanism by which the walls bounded by strings dissipate their 
energy since each string is being pulled in different directions by the N walls to which it is 
attached. 

The fraction r of wall energy that goes into non-relativistic or barely relativistic axions 
is at any rate severely restricted since r is much earlier than the time teg ~ 3 • 10^^ sec of 
equality between matter and radiation and, using Eqs. ( |6.3| ) and (|6.7D, 

r -10" — >10-M . 6.8 

\ vfia ) sec \ rUa J 

The fraction of the cosmological energy density in axions from domain wall decay 

Qt'"it)=rQ,.Ur)l^iy (6.9) 

should be less than about 1/2 when t = t^q. This implies the bound 

^ A ( "^n \ /secX-*^/^ ^ 4 / Vtlf, \'^/^ , , 

where we used Eq. ( |6.7|) to obtain the second inequality. 

However emission of gravitational waves is likely to be very effective because t^.^., given 
by Eq. (3.40), is not very much larger than r. Thus decay of the walls into gravitational 
radiation is in this case both necessary and likely. The third solution therefore predicts a 
peak in the present spectrum of gravitational waves at the frequency 

1 (^r\ ^ -.^-lOrr /Sec\l/2 



The energy density in these gravitational waves is of order 

/lO~^eV\ r 

P,.^{U) ~ 10-^7(^0) — , (6.12) 

\ rUa ) sec 

where P'y{to) is the energy density in the cosmic microwave background. 
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VII. CONCLUSIONS. 



We have discussed the fate of the walls bounded by strings which appear during the QCD 
phase transition in axion models with N — 1 assuming the axion field is not homogenized 
by inflation. We have argued that the main decay mechanism of these objects is radiation 
of barely relativistic axions and that the decay takes place during or soon after the QCD 

phase transition. 

We presented the results of our computer simulations of the motion and decay of walls 
bounded by strings. In the simulations, the walls decay immediately, i.e. in a time scale of 
order the light travel time. The computer simulations also provide an estimate of the average 
energy of the axions emitted in the decay of the walls: (oja) — Tnia when \/Xva/ma — 20. 

Because of restrictions on the available lattice sizes, the simulations are for values of 
Va/rria of order 100, whereas in axion models of interest Va/ma is of order 10^^. To address 
this shortcoming, we have investigated the dependence of {iOa) /ma upon \fXva/'ma and found 
that it increases approximately as the logarithm of this quantity. This is because the decay 
process occurs in two steps: 1) wall energy converts into moving string energy because the 
wall accelerates the string, and the energy spectrum hardens accordingly; 2) the strings 
annihilate into axions without qualitative change in the energy spectrum. If this behaviour 
persists all the way to \/~\Va/'ma ~ 10^^, then {ua)/ma ~ 60 for axion models of interest. 

To parameterize our ignorance of the decay process we introduced two parameters: the 
time ^3 when the decay effectively takes place and the average Lorcntz factor 7 = (tfa)/ma 
of the radiated axions at time ^3. We found that in first approximation the cosmological 
energy density of axions from wall decay does not depend upon ^3 whereas it is inversely 
proportional to 7. For 7 ~ 7 the contribution from wall decay is of the same order of 
magnitude as that from vacuum realignment. All three contributions - vacuum realignment, 
string decay and wall decay - were discussed in one unified treatment. 

We discussed the velocity dispersion of axions from wall decay (pop. II) and found that it 
is much larger, by a factor 10^ or so, than that of axions from vacuum realignment and from 
string decay (pop. I). This has interesting consequences for the formation and evolution of 
axion miniclusters. We showed that the QCD horizon scale density perturbations in pop. II 
axions get erased by free streaming before the time t^^ of equality between radiation and 
matter whereas those in pop. I axions do not. Pop. II axions form an unclustered component 
of the axion cosmological energy density which guarantees that the signal in direct searches 
for axion dark matter is on all the time. 

We speculated that the velocity dispersion of axions from wall decay may be large enough 
to be measured in a cavity detector of dark matter axions. This would provide a direct 
experimental handle on some of the issues raised by axion cosmology, e.g. the question 
whether inflation occurred below or above the PQ phase transition. 

Finally we discussed the cosmology of axion models with iV > 1 in which the axion 
domain wall problem is solved by postulating a small Upq(1) breaking interaction which 
slightly lowers one of the N vacua with respect to the others. There is very little room in 
parameter space for this to work but it is a logical possibility. We flnd that in this case the 
walls must decay into gravitational radiation of frequency of order 10~^°Hz. 
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FIG. 2. (a) The circle of axion field expectation values is divided into three equal parts with the 
CP conserving value a = in the middle of part 1. (b) Each cell of a hexagonal grid is randomly 
assigned 1,2 or 3. This results in a set of domain walls (thick lines), upgoing strings (open circles), 
and downgoing strings (filled circles), as described in the text. 



31 



r 




I — I 

horizon 
size 



FIG. 3. Same as in Fig. 2b except that the area covered is larger, the underlying grid has been 
removed and the walls have been smoothed. 
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FIG. 4. Speed of the string core as a function of time for the case l/rua = 400, y/X/nia = 10, 
and V = 0. 
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FIG. 5. Position of the string or anti-string core as a function of time for l/rria = 1000, 
A = 0.0002, D = 2896, and v = 0. The string and anti-string cores have opposite z. They go 
through each other and oscillate with decreasing amplitude. 
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(b) 

FIG. 6. Gradient (solid), kinetic (dotted), wall potential (long dash), string potential (short 
dash) and total (dot dash) energy as a function of time for the case l/rua = 1000, A = 0.0002, 
D = 2896, and v = (a), and for the case l/m„ = 500, A = 0.0032, D = 2096 and v = 0.25 (b). 
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(e) (f) 

FIG. 8. Energy spectrum at successive time intervals At = l/nia, with 

Iw = ^8log{w /nia) / ^og{wraax/'ma) + 2 and Wmax = v'8 + "^ai t'^^ "^^^e = 500, A = 0.0032, 
D = 2096 and v = 0.25. 
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FIG. 9. (w) as a function of time for l/m^ = 500, D = 2096, v = 0.25 and A = 0.0004 (solid), 
0.0008 (long dash), 0.0016 (short dash), 0.0032 (dot) and 0.0064 (dot dash). After the wall has 
decayed, (u) is the average energy of radiated axions. 
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FIG. 10. Same as in Fig. 9 except A = 0.0016, and v = 0.15 (sohd), 0.25 (long dash), 0.4 (short 
dash) and 0.6 (dot). 
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